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Abstract 

The viability of a possible cosmological scenario is investigated. 
The theoretical framework is the constrained next-to-minimal super- 
symmetric standard model (cNMSSM), with a gravitino playing the 
role of the lightest supersymmetric particle (LSP) and a neutralino 
acting as the next-to-lightest supersymmetric particle (NLSP). All the 
necessary constraints from colliders and cosmology have been taken 
into account. For gravitino we have considered the two usual produc- 
tion mechanisms, namely out-of equillibrium decay from the NLSP, 
and scattering processes from the thermal bath. The maximum al- 
lowed reheating temperature after inflation, as well as the maximum 
allowed gravitino mass are determined. 



1 Introduction 



There is accumulated evidence both from astrophysics and cosmology that 
about 1/4 of the energy budget of the universe consists of so called dark 
matter, namely a component which is non-relativistic and neither feels the 
electromagnetic nor the strong interaction. For a review on dark matter 
see e.g. |T]. Although the list of possible dark matter candidates is long, 
it is fair to say that the most popular dark matter candidate is the lightest 
supersymmetric particle (LSP) in supersymmetric models with R-parity con- 
servation [2]. The superpartners that have the right properties for playing 
the role of cold dark matter in the universe are the axino, the gravitino and 
the lightest neutralino. By far the most discussed case in the literature is the 
case of the neutralino (see the classic review [3j), probably because of the 
prospects of possible detection. 

However, the gravitino is another very interesting candidate for cold dark 
matter, since its interactions are completely determined by the supergravity 
lagrangian [4J, in contrast to what happens to the neutralino or the axino 
case, where the interactions depend on the chosen model. Unfortunately, 
gravitino belongs to the class of new exotic particles that can be potentially 
dangerous for cosmology, and it is therefore escorted by the so-called grav- 
itino problem [3] . The mass of the gravitino strongly depends on the SUSY- 
breaking scheme, and can range from eV scale to scales beyond the TeV 
region [6l [7J [8] . In particular, in gauge-mediated SUSY-breaking schemes [6] 
the gravitino mass is typically less than 100 MeV, while in gravity-mediated 
schemes [7J it is expected to be in the GeV to TeV range. Finally, it must 
be noted that there are hybrid models of gauge- and gravity-mediation, in 
which gravity provides sub-dominant and yet non- negligible contributions [S] . 
Therefore, according to the precise mechanism for supersymmetry breaking, 
the gravitino can be either stable or unstable, with the corresponding grav- 
itino cosmology. In general, the gravitino problem requires that the reheating 
temperature after inflation should be lower than 10 6 — 10 7 GeV [TQl [IT] , which 
poses serious difficulties to the thermal leptogenesis scenario [T2] . 

In the present work we want to consider gravitino dark matter in the con- 
strained next-to-minimal supersymmetric standard model (cNMSSM), as- 
suming that the lightest neutralino is the next-to-lightest supersymmetric 
particle (NLSP), and taking into account the two usual gravitino production 
mechanisms to be discussed later on, namely the out-of-equillibrium decays 
of the NLSP, as well as scattering processes from the thermal bath. 
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This article is organized as follows. In the next section we present the 
theoretical framework. In section 3 we discuss all the relevant constraints 
from colliders and from cosmology, and we show our results. Finally, we 
conclude. 

2 Theoretical framework 

In the present article we work in the framework of the constrained next-to- 
minimal supersymmetric standard model (cNMSSM). We assume that the 
gravitino is the LSP, while the lightest neutralino is the NLSP. The gravitino 
is stable and plays the role of cold dark matter in the universe, while the 
neutralino is unstable and it decays to gravitino. 

In what folows we review in short the particle physics model, namely the 
cNMSSM, as well as the gravitino production mechanisms. 



The most straightforward extension of standard model (SM) of particle physics 
based on SUSY is the minimal supersymmetric standard model (MSSM) [13J. 
It is a supersymmetric gauge theory based on the SM gauge group with the 
usual representations (singlets, doublets, triplets) and on Af = 1 SUSY. Ex- 
cluding gravity, the massless representations of the SUSY algebra are a chiral 
and a vector supermultiplet. The gauge bosons and the gauginos are mem- 
bers of the vector supermultiplet, while the matter fields (quarks, leptons, 
Higgs) and their superpartners are members of the chiral supermultiplet. 
The Higgs sector in the MSSM is enhanced compared to the SM case. There 
are now two Higgs doublets, H u ,Hd, (or Hi,H 2 ) for anomaly cancelation 
requirements and for giving masses to both up and down quarks. After elec- 
troweak symmetry breaking we are left with five physical Higgs bosons, two 
charged H ± and three neutral A, H, h (h being the lightest). Since we have 
not seen any superpartners yet, SUSY has to be broken. In MSSM, SUSY is 
softly broken by adding to the Lagrangian terms of the form 

• Mass terms for the gauginos Mi, M 2 , M 3 



2.1 Basics of cNMSSM 



M~g~9 



(1) 



• Mass terms for sfermions / 



(2) 
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• Masses and bilinear terms for the Higgs bosons H u , 

m\H\H u + m 2 H H ] d H d + Bfi{H u H d + h.c.) (3) 

• Trilinear couplings between sfermions and Higgs bosons 

AYf\Hf 2 (4) 

In the unconstrained MSSM there is a huge number of unknown parame- 
ters [H] and thus little predictive power. However, motivated by the grand 
unification idea, the constrained MSSM (CMSSM) assumes that gaugino 
masses, scalar masses and trilinear couplings have (separately) a common, 
universal, value at the GUT scale, like the gauge coupling constants do. 
CMSSM is therefore a framework with a small controllable number of pa- 
rameters, and thus with much more predictive power. In the CMSSM there 
are four parameters, mo,mu2, A ,tan(3, which are explained below, plus the 
sign of the \x parameter from the Higgs sector. The magnitude of //, as well 
as the B parameter mentioned above, are determined by the requirement for 
a proper electroweak symmetry breaking. However, the sign of /i remains 
undetermined. The other four parameters of the CMSSM are related to 

• Universal gaugino masses 

M^Mgut) = M 2 (M GUT ) = M 3 (M GUT ) = m 1/2 (5) 

• Universal scalar masses 

mj i (M GUT )=m (6) 

• Universal trilinear couplings 

A%(M GUT ) = A^Mgut) = A^Mqut) = A 6 tJ (7) 

tan(3 = — 

where v% , v 2 are the vevs of the Higgs doublets and Mqut ~ 
is the Grand Unification scale. 



(8) 

10 16 GeV 
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Unfortunately, the CMSSM suffers from the so-called fi problem [To] . 
This problem is elegantly solved in the framework of the next-to-minimal 
supersymmetric standard model (NMSSM) [16]. In addition to the MSSM 
Yukawa couplings for quarks and leptons, the NMSSM superpotential con- 
tains two additional terms involving the Higgs doublet superfields, Hi and 
H 2 , and the new superfield S, a singlet under the SM gauge group SU(3) C x 



where we take Hf = (H®,H[), H 2 = (H^H®), i,j are SU(2) indices, and 
ei2 = 1. In this model, the usual MSSM bilinear /i term is absent from 
the superpotential, and only dimensionless trilinear couplings are present in 
W. However, when the scalar component of S acquires a VEV, an effective 
interaction \iHiH 2 is generated, with /x = X(S). 

Finally, the soft SUSY breaking terms are given by [T7] 



Clearly, the NMSSM is very similar to the MSSM. Despite the similar- 
ities between the two particle physics models, the Higgs sector as well as 
the neutralino mass matrix and mass eigenstates in the NMSSM are more 
complicated compared to the corresponding ones in the MSSM. 

In particular, in the Higgs sector we have now two CP-odd neutral, and 
three CP-even neutral Higgses. We make the assumption that there is no 
CP-violation in the Higgs sector at tree level, and neglecting loop level effects 
the CP-even and CP-odd states do not mix. We are not interested in the 
CP-odd states, while the CP-even Higgs interaction and physical eigenstates 
are related by the transformation 



SU(2) L X U(l)y [I?] 





(10) 



h° - S ab H° 



(11) 
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where S is the unitary matrix that diagonalises the CP-even symmetric mass 
matrix, a,b = 1,2,3, and the physical eigenstates are ordered as m h o < 



m h o < m h o. 



In the neutrahno sector the situation is again more involved, since the 
fermionic component of S mixes with the neutral Higgsinos, giving rise to 
a fifth neutrahno state. In the weak interaction basis defined by ty° T = 

(^B° = = — 2A3, H®, H®, Sj , the neutrahno mass terms in the La- 



grangian are [T7] 



^ ass = -^(^T-Mx^ + H.C, (12) 



with Ai^o a 5 x 5 matrix, 

/ Mi -M z swO w cosP M z suxO W f&a.l3 

M 2 M z cos 9 W cos j3 —M z cos 9 W sin j3 

—Mz sin Ow cos (3 M z cos 8w cos (3 —As — Xv 2 

MzsinOwsmP — M z cos 9w sin/3 —As — \v\ 

\ -\v 2 -Xvi 2ks 

(13) 

The above matrix can be diagonalised by means of a unitary matrix iV 

^M^oN- 1 = diag(m^o, m^o, m^o, m^o, m^o) (14) 

where m^o is the lightest neutrahno mass. Under the above assumptions, the 
lightest neutrahno can be expressed as the combination 

Xi = N n B° + N 12 W° + N 13 H° + N 14 H° + N 15 S (15) 

In the following, neutralinos with Nf-^ > 0.9, or Nf 5 > 0.9, will be referred 
to as bino- or singlino-like, respectively. 

Similarly to the CMSSM, in the constrained next-to-minimal supersym- 
metric standard model the universality of mo, Ao,my2 at the GUT scale is 
again assumed, with the only parameters now being [18J 

tan/3, m Q , A ,m 1/2 , A, A k 
and the sign of the \x parameter can be chosen at will. 

We end the discussion on the particle physics model here, by making a fi- 
nal remark regarding the differences between the CMSSM and the cNMSSM. 
In the CMSSM the lightest neutrahno is mainly a bino in most of the param- 
eter space, and low values of m are disfavored because they lead to charged 
sleptons that are lighter than the neutrahno x?> while in the cNMSSM the 
lightest neutrahno is mainly a singlino in large regions of the parameter space, 



5 



thanks to which the charged LSP problem can be avoided [TS]. Furthermore, 
in the cNMSSM there are more mechanisms that contribute to the neutralino 
relic density [18]. 

2.2 Gravitino production 

For the gravitino abundance we take the relevant production mechanisms 
into account and impose the cold dark matter constraint [19] 

0.1097 < tt cdm h 2 = tt 3/2 h 2 < 0.1165 (16) 

At this point it is convenient to define the gravitino yield, Y3/2 = 
where n 3 / 2 is the gravitino number density, s = /i e //(T)^T 3 is the entropy 
density for a relativistic thermal bath, and h e jf counts the relativistic degrees 
of freedom. The gravitino abundance ^3/2 in terms of the gravitino yield is 
given by 

a** = ma8 ™ >V ' 2 = 2.75 x 10° r 3/2 (T 0) (17) 

where we have used the values 

T = 2.73K = 2.35 x 10~ 13 GeV (18) 
h eff (T ) = 3.91 (19) 
p cr /h 2 = 8.1 x 10~ 47 GeV 4 (20) 

The total gravitino yield has two contributions, namely one from the thermal 
bath, and one from the out-of-equillibrium NLSP decay. 

Y3/2 = Y3/2 + Y3/2 P ( 21 ) 

The contribution from the thermal production has been computed in [2U1 
I2"T| [22]. In [20] the gravitino production was computed in leading order in 
the gauge coupling g 3 , in [21] the thermal rate was computed in leading order 
in all Standard Model gauge couplings gY,92,93, and in [22] new effects were 
taken into account, namely: a) gravitino production via gluon — > gluino 
+ gravitino and other decays, apart from the previously considered 2—7-2 
gauge scatterings, b) the effect of the top Yukawa coupling, and c) a proper 
treatment of the reheating process. Here we shall use the result of [20] since 
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the corrections of [2H 122] do not alter our conclusions. Therefore the thermal 
gravitino production is given by 



- 29 x 10 " 2 (re) f 1 + 3) (22) 

or, approximately for a light gravitino, <C 

^-°- 13 (ira)fe) 2 (23) 

with tjIq the gravitino mass and m~ g the gluino mass. At this point it must 
be noted that all the relevant particles here (gluons, quarks, gluinos, squarks 
and gravitinos) are supposed to be in thermal equillibrium, and thus the 
reheating temperature after inflation should be at least 1 TeV. This is going 
to be important later on. 

The second contribution to the gravitino abundance comes from the decay 
of the NLSP 

n NLSP h 2 = QNLsph 2 (24) 

1 m NLSP 

with ttinlsp the mass of the NLSP, and flNLSph 2 the abundance the NLSP 
would have, had it not decayed into the gravitino. 



3 Constraints and results 

- Spectrum and collider constraints: We have used NMSSMTools (23], a 
computer software that computes the masses of the Higgses and the super- 
partners, the couplings, and the relic density of the neutralino, for a given set 
of the free parameters. We have performed a random scan in the whole pa- 
rameter space (with fixed /x > motivated by the muon anomalous magnetic 
moment), and we have selected only those points that satisfy i) theoretical re- 
quirements, such as neutralino LSP, correct electroweak symmetry breaking, 
absence of tachyonic masses etc., and ii) LEP bounds on the Higgs mass, 
collider bounds on SUSY particle masses, and experimental data from B- 
physics [2U [2S]. For all these "good" points the lightest neutralino is either 
a bino or a singlino, and contrary to the case where neutralino is the dark 
matter particle, here we do not require that the neutralino relic density falls 
within the allowed WMAP range mentioned before. 
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- As we have already mentioned, the total gravitino abundance, and not 
the neutralino one, should satisfy the cold dark matter constraint [TP] 

0.1097 < tt cdrn h 2 = tt 3/2 h 2 < 0.1165 (25) 

that relates the reheating temperature after inflation to the gravitino mass 
as follows 

0.11 = A(m^m- g )T R + ^^n NLS ph 2 (26) 

TTlNLSP 

For a given point in the cNMSSM parameter space, the complete spectrum 
and couplings have been computed, and we are left with two more free pa- 
rameters, namely the gravitino mass and the reheating temperature after in- 
flation. The gravitino mass is directly related to the SUSY-breaking scheme, 
while the precise range of values of the reheating temperature is crucial for 
the baryon asymmetry generation mechanism. The thermal production con- 
tribution cannot be larger than the total gravitino abundance, and for this 
we can already obtain an upper bound on the reheating temperature 

T fi <4.1xl0 9 f (™Y GeV (27) 

K ~ V100 GeV J \m~g ) K J 

Assuming a gluino mass rrig ~ 1 TeV, we can see that for a heavy gravitino, 
tuq ~ 100 GeV, it is possible to obtain a reheating temperature large enough 
for thermal leptogenesis. However, for a light gravitino, tuq ~ 1 GeV, we 
cannot obtain a reheating temperature larger than Tr ~ 10 7 GeV. 

- In scenarios in which gravitino is assumed to be the LSP, the NSLP 
is unstable with a lifetime that is typically larger than BBN time tsBN ~ 
1 sec. Energetic particles produced by the NLSP decay may dissociate the 
background nuclei and significantly affect the primordial abundances of light 
elements. If such processes occur with sizable rates, the predictions of the 
standard BBN scenario would be altered and the success of the primordial 
nucleosynthesis would be spoiled. BBN constraints on cosmological scenarios 
with exotic long-lived particles predicted by physics beyond the Standard 
Model have been studied [TDl EH]- Previous investigations have shown that 
the neutralino NLSP scenario with a gravitino mass m.3/2 > 100 MeV and 
a neutralino lifetime in the range (10 4 — 10 8 ) sec is already disfavored [27] . 
while the stau NLSP is still a viable scenario. The neutralino NLSP scenario 
can still be rescued if we avoid the stringent BBN constraints, namely if 
the neutralino lifetime becomes either larger than the age of the universe or 
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lower than the BBN time. The first possibility is realized in the degenerate 
gravitino scenario [2S], where the neutralino is extremely long-lived, and the 
only constraint comes from the cold dark matter bound. 

For neutralino NLSPs, the dominant decay mode is \ ~ ^ 1^, for which 
the decay width is [27J [29] 



| N n cos 9 W + N12 sin W | 2 m 
48vrM2 m 



rn 



a 



mi 



1 + 3 



m?. 



(28) 



where M* is the Planck mass, m x is the neutralino mass, and 6\y is the weak 
angle. This decay contributes only to EM energy. If kinematically allowed, 
the neutralino will also decay to gravitino and Z boson, or gravitino and light 
standard model Higgs boson h. The leading contribution to hadronic energy 
is from x ~> ZG, hG. These decays produce EM energy for all possible Z 
and h decay modes (except Z — > vv\ but they may also produce hadronic 
energy when followed by Z, h — > qq. The decay width to Z bosons is J27J [29] 



r( x -> ZG) 



— Nu sin 9w + N12 cos 9 



w 



m" 



1 - 




m 



-F(m x ,m d ,m z ) 



G 



1 + 3- 



m, 



m 



mi 



rn 



G(m x ,m^,m z ) 



29) 



where 

F(m x ,m 6 ,m z ) 



m^ + m z 



rn 



x 




m. 



1/2 

(30) 



G(m x ,m d ,m z ) = 3 + (-12 + ^+ ^f-^f(3-^f) (31 ) 

m 6 x \ m x J m x m x \ m l x ' 



with m z — 91 GeV the mass of the Z boson. 

The decay width to the Higgs boson is [271 EH] 

w \N 13 S U + N U S 12 + Ni 5 S 13 \ 2 m 5 
T(x^hG) = to-.*,-?. —F{m x ,m d ,m h ) 



48ttM 2 

2 



m- 



— \ 1 + — - -%H(m x ,m d , m h ) 
m xJ \ m xJ m x 



(32) 
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where are three of the components of the mixing matrix in the Higgs 
sector for the CP-even mass eigenstates, F is as given in (I30p . and 



rrin m % . m % m % mi 



H{m x ,m d ,m h ) = 3 + 4^ + 2-f + 4-f + 3-f + -} 

m y mi m% m% m% 

A A A A 

-4f3 + 2^ + 34) (33) 

with nth the mass of the Higgs boson. Therefore, the neutralino lifetime is 
given by 

t = f (34) 
r = T(x^lG) + T(x^ZG) + T( X ^hG) (35) 

Given these two-body decay widths, the resulting values for the energy release 
parameters are 

Bl M * 1 (36) 

2 2 

m x — m^ 

£ ™ = 2m x (3?) 

BX „ r( x -> ZG)Bj d + T( X -» ^G)< d + r( x -> ggG) (3g) 
had r( x ^ 7 G) + r( x ^zG) + r( x ^/iG) 

£had ~ 2m x ' ldJj 

where -B^ ad ~ 0.9, B^ ad m 0.7, and the three-body decay T(x — > <?<?G0 ~ 
io- 3 r [2Z1. 

- Finally, it must be noted that for long neutralino lifetimes, t^lsp > 
10 7 sec, in addition to BBN constraints there are strong bounds from the 
shape of the cosmic microwave background (CMB) black-body spectrum [281 
130] . However, in our investigation we have found that the neutralino lifetime 
is always r NLS p < 10 2 sec, for gravitino masses m^ < 1 GeV, and therefore 
we do not need to worry about these bounds from the CMB shape. 

Our main results are summarized in the figures below. Before starting to 
discuss the figures, let us first make a few comments. The precise neutralino 
composition depends on the values of the coefficients Nu in ( Fl5|) . which in 
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turn depend on the values of the free parameters of the model. Roughly, for 
large coupling A = 0.1 — 0.5 the neutralino is mainly a bino, while for small 
coupling A< 1, the neutralino is mostly a singlino. We refer the interested 
reader to e.g. (18] for the relevant discussion. Furthermore, the neutralino 
lifetime is determined by the three decay channels to gravitino plus photon 
or Z boson or Higgs boson. These partial decay rates depend on the available 
phase space (masses) as well as the couplings (composition coefficients Nu). 
Thus in the bino case, in which Nu ~ 1 and the rest of the coefficients 
are very small, the decay rate to gravitino and Higgs boson is negligible, 
while in the singlino case, in which N15 ~ 1 and the rest of the coefficients 
are very small, the decay rate to gravitino and Z boson is negligible. The 
decay channel to gravitino and photon gives the main contribution, while the 
decay rate to gravitino and Z boson (in the bino case) or to gravitino and 
Higgs boson (in the singlino case) modify the neutralino lifetime by a factor of 
twenty or fifty per cent respectively. Finally, for a given point in the cNMSSM 
parameter space, the neutralino lifetime is a function of the gravitino mass 
only. Imposing the BBN constraints we find the maximum allowed gravitino 
mass, and from the cold dark matter bound we can determine the maximum 
allowed reheating temperature. 

We can now turn to the figures where we discuss the two cases (bino 
or singlino) separately. The first three figures correspond to the bino case, 
while the last two figures correspond to the singlino case. For the bino case, 
it is important to notice that the neutralino relic density can take values 
larger than the usual ones by two orders of magnitude. The reheating tem- 
perature decreases with the neutralino relic density, and takes larger values 
for very low neutralino relic density. Figure 3 shows the maximum allowed 
reheating temperature after inflation versus the maximum allowed gravitino 
mass, both in GeV. Although it cannot be seen directly from the figures, the 
maximum possible gravitino mass in the bino case is rriQ ~ 1 GeV, and the 
corresponding reheating temperature is Tr ~ 10 7 GeV. Therefore, we see 
that a) the gravitino in this scenario must be much lighter than the rest of 
superpatners, and b) the reheating temperature after inflation is not large 
enough for thermal leptogenesis. We remark in passing that gravitino masses 
of the order of 1 GeV can be obtained in hybrid supersymmetric models of 
gauge- and gravity-mediation, in which gravity provides sub-dominant and 
yet non-negligible contributions [9]. 

For the singlino case, we show in figure 4, the gravitino mass in GeV ver- 
sus neutralino relic density, and in figure 5 the maximum allowed reheating 
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temperature versus gravitino mass, both in GeV. This time the neutralino 
relic density is even larger than before, and gravitino now must be extremely 
light. This is due to the smallness of the coefficients Nn,Nu in the decay 
rate to gravitino and photon. For the same lifetime as before, the gravitino 
mass must be several orders of magnitude lower than in the bino case. The 
last figure shows that in the singlino case the reheating temperature can- 
not be larger than about 200 GeV. However, this value is much lower than 
the minimum value required for the computation of the gravitino thermal 
production, and therefore we conclude that this scenario must be excluded. 
We can understand these features as follows. First, recall that the WMAP 
bound for the cold dark matter abundance relates the gravitino mass to the 
reheating temperature, and we can obtain an upper bound on the reheat- 
ing temperature for a given gravitino mass. From equation (27), and for a 
gluino mass m~ g ~ 1 TeV, we see that when the gravitino is extremely light, 
rriQ ~ 10 -6 GeV, the upper bound on the reheating temperature becomes 
T R ~ 41 GeV. We then need to understand why the gravitino becomes so 
light in the singlino case. Let us assume that we have in the parameter space 
a point that corresponds to the bino case, another point that corresponds to 
the singlino case, and that the Higgs mass, superpartner masses, as well as 
the neutralino abundance are the same for the two points. The only thing 
that is different is the composition coefficients for the neutralino. In the 
bino case the first coefficient is practically unity and the rest tiny, while in 
the singlino case the last coefficient is almost unity and the rest negligible. 
The BBN constraints determine the maximum possible gravitino lifetime, 
which is given essentially by the photon decay channel. We thus have for the 
neutralino lifetime 

m% 

T ~W^ m 

Therefore if in the singlino case |-/V n | ~ 10~ 6 , the gravitino mass becomes as 
low as mg ~ 10~ 6 GeV. 



4 Conclusions 

In the framework of the cNMSSM, we have considered a possible cosmological 
scenario with the gravitino LSP and the neutralino NLSP. The gravitino 
is stable and plays the role of cold dark matter in the universe, while the 
neutralino is unstable and decays to gravitino. We have taken into account 
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the relevant gravitino production mechanisms, which are i) the NLSP decay, 
and ii) scattering processes from the thermal bath. Our results can be seen 
in the figures. We have found that i) the gravitino is necessarily very light, 
and ii) the reheating temperature after inflation is two orders of magnitude 
lower than the temperature required for thermal leptogenesis. The singlino 
scenario must be excluded, while in the bino case it is possible to have a 
gravitino in the gravity-mediated SUSY-breaking scheme. 
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Figure 1: Gravitino non-thermal production versus neutralino relic density 
for the bino case. 
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Figure 2: Reheating temperature (in GeV) versus neutralino relic density for 
the bino case. 
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Figure 3: Reheating temperature versus gravitino mass (both in GeV) for 
the bino case. 




Figure 4: Gravitino mass (in GeV) versus neutralino relic density for the 
singlino case. 



18 




Figure 5: Reheating temperature versus gravitino mass (both in GeV) for 
the singlino case. 
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